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ABSTRACT 



Context. The quiet Sun presents magnetized plasma whose field strengths vary from zero to about 2 kG. The probability density 
function of the magnetic field strength B efficaciously describes the statistical properties of the quiet Sun magnetic field. 
Aims. We simulate the dynamics and the evolution of quiet Sun magnetic elements to produce a probability density function of the 
field strengths associated with such elements. 

Methods. The dynamics of the magnetic field are simulated by means of a numerical model in which magnetic elements are driven 
passively by an advection field characterized by spatio-temporal correlations that mimick the granulation and mesogranulation scales 
observed on the solar surface. The field strength can increase due to an amplification process that occurs where the magnetic elements 
converge. Starting from a <5-like probability density function centered on 6 = 30 G, we obtain magnetic field strengths of up to 2 kG 
(in absolute value). To derive the statistical properties of the magnetic elements, several simulation runs are performed. 
Results. Our model is able to produce kG magnetic fields in a time interval of the order of the granulation timescale. The mean 
unsigned flux density and the mean magnetic energy density of the synthetic quiet Sun reach values of (6) ^ 100 G and (6")''- ^ 
350 G respectively in the stationary regime. The derived probability density function of the magnetic field strength decreases rapidly 
from 6 = 30 G to 6 ~ 100 G and has a secondary maximum at 6 = 2 kG. From this result, it follows that magnetic fields > 700 G 
dominate the unsigned flux density and magnetic energy density, although the probability density function of the field strength reaches 
a maximum at B ~ 10 G. 

Conclusions. A photospheric advection field with spatio-temporal correlations, driving the magnetic elements, and reduced magnetic 
amplification rules are able to create a realistic probability density function of the quiet Sun magnetic field. It has been found that they 
naturally produce an excess of magnetic fields around 2 kG if an upper limit is imposed on the field strength. 
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, When observing the Sun in full disk magnetograms, magnetic 

■ concentrations define the network pattern associated with 
the supergranular motions of photospheric plasma, while the 

, interior of the network pattern appears to be unmagnetized. 

' We know that weak polarization signals are measur- 
able almost everywhere in the interior o f network cells 
([Livingston & HarvevI 119751; ISmithsonI Il975h . These signals 
are associated with the internetwork component of the solar 
magnetic field, also called the quiet Sun magnetic field, which 
covers more t han 90% of the solar surface, independent of the 
solar activity dHarvev- An gld [T993h . Due to its large diffusion 
onto the solar surface, the quiet Sun magnetic field may be 
useful for understanding the physics of the photospheric mag- 
netism. In fact, it could hold a large fraction of the unsigned flux 
dens ity and mag netic en ergy density of the solar photosphere 
(e.g. 'Unno"1959'; 'Stenfld ll982l: ISanchez Almeidalll998i 120041: 
[Schrijver & Title 2003). 

The weak polarization signals measured in quiet Sun regions 
can be interpreted as the result of the linear combination 
of the polarization emerging from discrete magnetic flux 
tubes, which are typically smaller than the angular resolution, 
that fill the solar atmosphere with a complex to pology (e.g. 
lEmonet & Cattaiiedll2001al : ISanchez Almeida et al. ||2003). 
Knowledge of the quiet Sun has been enhanced sig- 
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2002; Sanchez Almeida 2004; Trujillo Buen o et al 
2004; Manso Sainz et al. 2004|; ISanchez Almeidal l2005l 
.Martmez Gon zalez et al., i2006bl). a nd numerical simula- 
tions (e.g. Cattaneo 119991; lEm onet & Cattaneo 2001 
Stein & Nordlund 2002; Vogler 2003; VogleretalJ 120051 
Stein & Nordlund 2006 ; Vogl er & S chussler 2 007) . 
These improvements were crucial in characterizing the quiet 
Sun magnetic field and its different observational parameters. 
One parameter, the fraction of solar surface occupied by a 
certain field strength, is a very powerful tool. In fact, this can be 
trivially translated into the probability density function (PDF) of 
the magnetic field strength, P(B), which is useful in describing 
the statistical properties of a quiet Sun region. For example, 
calculating the first and the second momentum of P{B) provides 
information about the mean unsigned flux density and the mean 
magnetic energy density (apart from a constant factor l/B^n), 
respectively. Using a PDF assumes implicitly a continuum for 
values of the magnetic field strength. The limit imposed by 
the pressure equilibrium between the magnetic flux tubes and 
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photospheric plasma also defines the variability range for quiet 
Sun magnetic fields to be between G to about 2 kG. 
[Lin (1995) computed the PDF in active and quie t Sun r egions 
using Fe I /II 5648 and Fe 1^15652 lines, while iLitesI (l2002h 
obtained it from the analysis of network regions and their 
surroundings using Fe i lines at 6300 A. Both authors found a 
bimodal distribution with a dominant contrib ution at ^ 150 G 
and a secondary peak at ^ 1.5 kG. In contrast, iColladosI (1200 ll) . 
studying q uiet Sun regions using the same infrared Unes as 
iLinI (1199 5'), found an exponential decline in the magnetic flux 
density histogram from a pe ak at ^ 200 G and no counts for 
magnetic flux density ^ 1 kG. lDommguez Cerdena et"al ] (l2006h 
developed a set of PDFs for the quiet Sun magnetic fields 
that was consistent with observations and magnetoconvection 
simulations by using information from both Hanle and Zeeman 
observations; their study created substantial interest because 
the functions increased for strong fields, just before a cut-off' 
at B ^ 1.8 kG. This characteristic implies that strong fields are 
important because, even when present in a small fraction of 
the photosphere, they dominate the unsigned flux density and 
magnetic energy density. On the othe r hand, the magnetic field 
dist ributions proposed by Martinez Gonzalez et al.l (l2006albl) 
and iMartmez Gonzalez et al., (2008.) , which were derived from 
simultaneous infrared and visible observations of Fei lines, 
do not present any rise for strong fields. The inversion of the 
FeU630 1 and Fe 1/16302 profiles observed by the HINODE 
satellite (lOrozco Suarez et aLll2007l) confirms this result. 
Magnetoconvection simulations provide a wide variety of 
quiet Sun magnetic field PDFs to be compared with those 
derived from observations. Simulations show, for exam- 
ple, exp onential decay s (Cattaneo 19991 IStei n et al] 120031; 



IStein & N ordlund 120061; JVogler & Schussler bim odal 
behaviors tVoder et al.( l2003i: iVogler & Schusslej l2003l) . or 



1 kG 



an "exponenti al-hke core and a modest peak" for B 
dSteiner 2003). 

An important step is to ascertain the presence or absence of a 
bump for magnetic fields ~ 1 kG in quiet Sun magnetic field 
strength PDFs; since the magnetic flux and magnetic energy 
both scale as powers of the field strength, kG magnetic fields 
should play a crucial role in defining the global properties 
of quiet Sun magnetic fields. The connectivity between the 
lower and the upper solar atmosphere may indeed be defined 
by t he larger field strength concentrations in the photosphere 
([Dommg uez Cerdena etal. 2006). 

ISanchez Almeidal ( f2007h indicated that a bump in the PDF 
for quiet Sun fields ~ kG is to be expected if "a magnetic 
amplification mechanism operates in the quiet Sun". The author 
derived an equation for the temporal evolution of the PDF, 
whose shape can be modified by magnetic amplification and 
photospheric convection. In the stationary state, the model is 
characterized by a P{E) with a bump for strong fields. 
In our model, we implement the dynamics of a high number 
of photospheric magnetic elements, driven by a photospheric 
velocity field in which a large-scale organization (mesogranu- 
lation) is naturally derived from the cooperation of small-scale 
flows (granulation). The model also takes into account simpli- 
fied processes that are able to amplify and limit the magnetic 
field strength. In § |2] we describe the model, in § |3] the results 
of the simulation are presented, while in § IH they are discussed. 
Finally, in §|5] the conclusions are outlined. 
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Fig. 1. Illustrative representation of \B\ in the computational do- 
main. Black areas represent G regions; blue and sky-blue rep- 
resent ~ 10 G regions; green represents ~ 100 G regions; yellow 
and red represent ~ kG magnetic concentrations. The position of 
concentrations with B ~ kG correspond to the position of stable 
downward plumes organized over mesogranular scales. 



2. The model 

The model simulates the dynamics and evolution of photo- 
spheric magnetic elements driven by a photospheric velocity 
field characterized by spatio-temporal correlations. To repro- 
duce these correlations, we simulated the phot ospheric vel oc- 
ity field in a way similar to that proposed by iRastI (l2003h ; it 
was computed by an n-body simulation that realizes the inter- 
action between downflow plumes by considering only their mu- 
tual horizontal interaction, as in Rast ( 1998). The movement of 
each plume was determined by the advection produced by the 
other plumes. This model spontaneously creates stable down- 
ward plumes, organized on a mesogranular pattern, from the in- 
teraction of granular scale flows. The ratios between the meso- 
granular and granular spatio-temporal scales, over which the 
model self-or ganizes, agree with th e ratios measured in the solar 
photosphere. Viticchie et al.l (l2006l) utilized such a velocity field 
to drive magnetic loop footpoints in a reconnection model that 
reproduces the observed statistical properties nanoflares. 

In the present work, the simulations start with the advection 
field only, without any magnetic element in the domain. When 
the velocity field has reached a stationary behavior, i.e. after the 
onset of mesogranular scale flows, the time t is set to zero and 
No - 500 magnetic elements are distributed at random positions 
over the computational domain. The magnetic elements have an 
associated field strength B,„ = 30 G with random orientation 
(+1) and size equal to the spatial resolution element. Hereafter, 
B represents the absolute value of the magnetic field strength. 
We choose such a value for Z?,„ since it is approximately equal to 
the mini mum magnetic flux density measurable nowadays (e.g. 
Truiillo Bueno et al.l2004l ; lBommier et al .120051; ICarroU & Kovi 
12008 ). We allow the magnetic elements to be transported pas- 
sively by the photospheric velocity field towards different plume 
sites. At each plume site, a certain number of magnetic elements 
converge under the action of the advection field. When two mag- 
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netic elements, associated with Bi and B2, overlap in the same 
resolution element, they can produce the following two results: 

1. If the two elements have opposite orientations, they are 
substituted by a new magnetic element at the same site with an 
associated B = \Bi - B2I and the orientation of the strongest ele- 
ment. The total flux is conserved such that the spatial size of the 
new element remains as one resolution element. In the particular 
case of Bi = B2, the elements cancel the effect of each other 

2. If the two elements have the same orientation, they are 
substituted by a new magnetic element in the same site with an 
associated B = B\ + B2 and the same orientation. The total flux 
is again conserved and the spatial size of the new element is one 
resolution element. 

This process implicates an amplification of the magnetic field 
strength. Different mechanisms have been proposed to concen- 
trate magnetic fields in quiet Sun regi ons u p to the kG range 
(e.g. I WeissI [19661: iPark er 1978'; Spruit 1979: Sanc hez Almeidal 
l200lh . As reported by Sanchez Almeida (2007) , all mecha- 
nisms are abruptly diminished at the maximum field strength 
imposed by the gas pressure of the quiet Sun photosphere. 
Therefore, we place the constraint that the amplification pro- 
cess cannot create field strengths greater than B/„„ = 2 kG. 
We can interpret B;„„ as a representative value of th e max- 
imum field 



strengths for quiet Sun regi ons ( e.g. iParkei 



1978; Spruit 1979; Grossmann-Doerth et al. 1998; Stein et al 
2003; Domfnguez Cerdeiia et al. 2006; Sanchez Almeida 2007). 
Whenever the amplification process would produce B > Bii,„, 
the magnetic elements do not interact and remain separated. 
We consider implicitly that the amplification mechanism op- 
erates on the same timescales as the granulation timescale 
(~ 10 min). This hypothesis is consiste nt with the 
amplification timescale s repo rted by iParkeJ ( Il978l) and 
iGrossmann-Doerth et al.l ( Il998l) . Up to the equipartition value 
of about 500 G, the amplification can be ascribed to the granular 
flows and five minutes are sufficient to produce kG fields from a 
uniform distribu tion of 400 G. 

As proposed bv lSteinerl ( l2003h . starting from an initial PDF of 
Pt=o(B) = 6(B - Bin), the amplification process will broaden the 
PDF and extend it to B/„„. 

To simulate the emergence of new magnetic elements on the so- 
lar surface, we continuously add randomly located elements of 
random orientation and B,„ magnetic field strength. This injec- 
tion process maintains the total number of magnetic elements 
~ A^o and compensates for the reduction in the total number due 
to the concentration process. 

Figure[T] shows an illustrative representation of a snapshot of the 
simulation in the stationary regime: we represent the \B\ pattern 
emerging on the computational domain. 

3. PDF of the quiet Sun magnetic field strengths 

To derive a PDF of the quiet Sun magnetic field strengths, we 
performed ten simulation runs with different advection fields and 
different initial distributions of magnetic elements. Each simula- 
tion reproduced the evolution of the quiet Sun for an equiva- 
lent time of about 2600 min. We considered the behavior of the 
model by analyzing its average maximum field strength max{B), 
mean unsigned flux density {B}, and mean magnetic energy den- 
sity {B^)l%n, as functions of the mean evolution time f: 



< B >= BP(B)dB, 
Jo 

< B^ >= B^P{B)dB, 



2000 - 
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Fig. 2. Average maximum field strength versus evolution time t. 
The function is obtained by averaging the results from rt/mN = 
10 distinct simulation runs and the error bars are calculated to 



equal 



. The horizontal dashed line represents the maximum 



value allowed for the magnetic flux density (B/,,,,), while the ver- 
tical dot-dashed line indicates the instant at which the model pro- 
duces the first kG element. The plot shows only a part of the total 
time domain (f ^ 2600 min). 
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Fig. 3. Mean probability density function (P(B)) for the mag- 
netic field strength in the stationary regime. (P(B)) has been ob- 
tained averaging tijoT - 40 non-correlated PDFs and the er- 
ror bars are calculated to equal -==. The histogram bin-size is 

30 G. The dashed line represents the upper limit for the mag- 
netic amplification process B/,>„. The maximum at B = G is 
produced by the fraction of the domain free of magnetic ele- 
ments. 



in a similar way to'Dommg uez Cerdena et al] (l2006h . 
The model creates magnetic concentrations with B = 1 kG in 
about 30 min and B ^ B;„„ in about 200 min (Fig.|2]l. 
The quantity max{B) is unsuitable for describing the stationary 
regime of the magnetic element system because its value is lim- 
ited. As discussed in § |2] we also expect to obtain a PDF span- 
ning the entire domain B = - B/,„,; the stationary regime must 
therefore be identified referring to the integral properties of the 
PDF, such as the mean unsigned flux density and the mean mag- 
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Fig. 4. Mean probability density function {P(B±)) for the signed 
magnetic field strength in the stationary regime. The histogram 
bin-size is 30 G. The vertical dashed lines represent the upper 
and lower limits of the magnetic amplification process ±B/„„- 
The maximum at Z? = G is produced by the fraction of the 
domain free of magnetic elements. 

netic energy density. {B) and (B^y^^ saturate exponentially at 
(B)sat - 100 G and {B^)l'J - 350 G, respectively. We inter- 
pret this saturation as an indication that the entire system has 
reached the stationary regime. In this regime, we compute {PiB)) 
by averaging forty non-correlated PDFs. The computed {P(B)} 
is shown in Fig.|3] 

Alternatively, we could compute the mean probability density 
function (P(B±)) for the signed magnetic field strength by con- 
sidering the direction of the magnetic fields (shown in Fig.|4]i. As 
expected, the {P{B±)} is quite symmetric with respect to zero. 

4. Discussion 

In the previous section we presented the results of several simu- 
lations of our reduced model. We now focus on some particular 
points in detail. 

The maximum field strength, plotted in Fig. |2] as a function of 
t indicates that the timescale over which kG fields are produced 
is approximately 30 min. This is the same timescale over which 
th e mesogranular patter n emerged in the numerical simulation 
of lViticchie eTai](l2006l) . 

Differ ent descriptions of the magnetic am plific atio n pro- 
cess (iGrossmann-Doerth et all Il998t ISanchez Almeida! 1200 Ih 
achieved similar results: weak magnetic-field structures, in a 
downflow environment, de velop field strengths ~ kG within a 
few minutes. Berrilli et al.| (|2005| ) reported that the organization 
of structures on mesogranular boundaries is rapid with a typical 
timescale of ^ 10 min. 

The analysis of the position of strong field magnetic elements 
in our simulation revealed that these are superimposed onto sta- 
ble downward plumes on mesogranular scales (Fig. [T]i. Such a 
correlation between strong magnetic field s and the mesogran- 
ular ad vection scale is in agreement with iDommguez Cerdefia 
(I2003h . who found a correlation between mesogranular photo- 
spheric flows and the pattern observed in magnetograms. 
The symmetry of the {P{B+)}, reported in Fig. |4l confirms 
the equal probability of having concentrations of positive and 
negative direction. This originates in both the initial condi- 



tions and the new element addition criteria applied during 
the evolution. Similar symmetries have been found in simu- 
lations (e.g., i Vogler & Schiissler 2007J and in observations of 
iMartmez Gonzalez et al] (l2006ai |2008|) . In the latter works the 
authors pointed out that such a symmetry may imply a local dy- 
namo origin for the internetwork magnetic field, that is a mag- 
netic field controlled by solar granulation. 
The sharp count decay found at B = B/„„ is obviously due to 
the "interaction turn off" produced by (§ |2]i. Other PDFs 
reported in t he literature show a smooth decay of the peak for 
B ^ 1 .8 kG jDomfnguez Cerdefia et a n i2006l:ISanchez Almeidal 
l2007 b. as expected for natural PDFs. We could reproduce this 
smooth decay of {P(B)) for B > B;„„ by modifying the magnetic 
element interaction probability, but it is not within the scope of 
this work to reproduce these fine details of the PDF. 
To verify that the stationary regime has been reached, we com- 
pared the {P{B)) with four mean PDFs, P,{B), representative 
of four different instants of the simulation during the station- 
ary regime (Fig.|5]l. Each Pt{B) represents the mean of ten PDFs 
at a given time t. We consider the four Pt{B) to be uncorrected 
because they are separated by more than 140 minutes, a time 
interval that is longer than three times the mean mesogranular 
lifetime. As shown in Fig. |5] the deviation between the {P{B)) 
and the four Pt{B) shows no time dependence and the root mean 
square deviations of all four Pt(B) with respect to the {P{B)) are 
approximately one part in a thousand. 

The values <B),«, ^ 100 G and {B-)\[^, ^ 350 G, reached in 
the stationary regi me, are in good agreement with the results of 
Dommguez Cerde fia"et al.l (|2006). The percentages of unsigned 
flux density and magnetic energy density for B > 700 G are 
- 65% and ^ 97% respectively, indicating a dominant role 
for kG fields. Th e percentages that we found are higher than 
those reported by 'Dommg uez Cerdefia et al.l (l2006l) : these dif- 
ferences can be explained by referring to the steeper decrease 
in the {P{B)) for the - 700 G interval with respect to the log- 
normal trend found by these authors. 

The (P(B)) presents a maximum at B = B„, in agreement with 
iDomrnguez Cerdefia et al.l (l2006h . This result originates in the 
injection value for the magnetic field strength Z?„, = 30 G. 
At Z? ^ Biim, the PDF p r esents a secondary maximum, a s 
fou nd by [S tein et al.1 (l2003|), iDommguez Cerdefia et al] (|2006|) . 
and [Sanchez Almeida (20071^ 

In a similar way to iSanchez Almeidal ('2007'), we may conclude 
that this secondary maximum in the {P{B)} is due to the upper 
limit placed on the amplification process at B;„„. When a mag- 
netic amplification process operates and a limit is placed on this 
process, an accumulation of field strength about the limit is ex- 
pected. All assumptions of the model contribute in shaping the 
{P{B)}: the strong convergence regions of the granular advection 
pattern create suitable conditions for amplifying the magnetic 
field and, at the same time, support the amplification process by 
gathering new magnetic elements to counteract the cancellation 
process. The secondary maximum of the PDF is produced by 
both the spatial and temporal correlations of the velocity field, 
which advect the magnetic elements, and the upper limit of the 
amplification process at B/,„,. 

5. Conclusions 

We have presented a model that reproduces the dynamics and 
evolution of photospheric magnetic elements. These elements 
are passively driven by an horizontal velocity field character- 
ized by spatio-temporal correlations that agree with those ob- 
served in the solar photosphere. The description of the velocity 
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Fig. 5. Mean probability density function {P(B)) (solid line, 
same as FigO) compared with four Pi(B) (dotted lines) calcu- 
lated at different instants of the simulation. Each Pi{B) repre- 
sents the mean of ten PDFs from ten simulation runs at time 
t. Top leff.t - 2160 min, top right t - 2308 min, bottom left 
t - 2464 min and bottom right t - 2593 min. The histogram 
bin-size is 30 G. The vertical dashed lines represent the upper 
and lower limits of the magnetic amplification process ±B/,„,. 
The maximum at B = G is produced by the fraction of the 
domain free of magnetic elements. 



field follows|Ras3 (Il998ll2003h . In their evolution, the magnetic 
elements interact and increase their magnetic field strengths 
through an amplification process assumed to work on a timescale 
comparable with the granular timescale and able to recreate 
kG concentrations. The injection scale, constrained to be B,„ - 
30 G, and the upper limit at B/„„ = 2 kG a re both derived from 
observations and physical constraints (e.g. iTruiiUo Bueno et al] 
[2004; Parker 1978). 

We investigated the statistical properties of the field strengths as- 
sociated with the magnetic elements using a series of numerical 
simulations. The main conclusions are: 

1 . The model produces kG fields in a time interval of the or- 
der of the mesogranulation timescale. These strong fields are 
organized in a pattern superimposed on the mesogranulation 
pattern of the advection field. 

2. In the stationary regime, the PDF of the magnetic field 
strength has a stable shape and has a secondary maximum 
that corresponds to the upper limit set at 2 kG. 

3. The associated mean unsigned flux density and mean mag- 
netic energy density in the stationary regime are (B),^, - 
100 G and {B^)l'J ^ 350 G, respectively. 

4. The magnetic element system described by this PDF is 
strongly dominated by kG magnetic fields: the fractions of 
unsigned flux density and magnetic energy density for B > 
700 G are ^ 65% and ^ 97%, respectively. 

The quiet Sun is certainly more complex than the model pre- 
sented here. However, a more simplified approach, based on 
the advection field proposed by Rast (2003), reproduced the ob- 
served distributions for the r eleased magnetic ene rgy and waiting 
times of nanoflare events (i Viticchie et alj|2006i) . In the present 
work, we showed how the same approach can also be applied 
to reproduce some of the statistical properties of the quiet Sun 
magnetic field. 
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